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Abstract 

The relativistic theory of elasticity is reviewed within the spherically symmetric context 
with a view towards the modeling of star interiors possessing elastic properties such as the 
ones expected in neutron stars. Emphasis is placed on generality in the main sections of 
the paper, and the results are then applied to specific examples. Along the way, a few 
general results for spacetimes admitting isometrics are deduced, and their consequences 
are fully exploited in the case of spherical symmetry relating them next to the the case in 
which the material content of the spacetime is some elastic material. This paper extends 
and generalizes the pioneering work by Magli and Kijowski Magli [2] and [3], and 
complements, in a sense, that by Karlovini and Samuelsson in their interesting series of 
papers [3], [5] and [^. 



1 Introduction 



The interest of a relativistic theory of elasticity is twofold; on the one hand there is 
its purely theoretical interest, namely that of providing a relativistic extension of a 
well-known (and very fruitful) classical theory; on the other hand and on theoretical 
grounds, it is expected that neutron stars possess a solid crust with elastic properties 
which may help explaining certain observational issues (see |1] for a thorough account 
of these and other related features). Further, anisotropy in pressures is a phenomenon 
occurring in many situations of equilibrium which are of interest in astrophysics and 
whose corresponding dynamics has been thoroughly studied (see for instance [7] and 
references therein); the assumed point of view however has been an heuristic one, 
without providing mechanisms explaining how the anisotropy in pressures may arise 
and using instead ad hoc assumptions. Magli and Kijowski [1] and Magli [3] have 
shown that in spherical symmetry, the anisotropy in pressures arises quite naturally 
as the relativistic extension of the classical (non-relativistic) non-isotropic stress in 
elasticity theory. See also [8] for an excellent study of the static case in spherical 
symmetry, where existence theorems for regular solutions near the center are proven 
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under rather mild, physically meaningful, hypotheses. Beig and Schmidt [9] have 
shown that, in general, the field equations for elastic matter can be cast into a first- 
order symmetric hyperbolic system and that as a consequence, local-in-time existence 
and uniqueness theorems may be obtained under various circumstances. 

The aim of this paper is to extend and generalize the work presented in [1] and [3] 
as well as to set up a set of mathematical tools and equations that may facilitate 
the obtention of exact solutions to Einstein's Field Equations (EFE) describing the 
interior of elastic materials and satisfying the Dominant Energy Condition (DEC). 

The paper is organized as follows: in the next section we provide a brief account of 
the theory of relativistic elasticity, much along the lines followed in pp, [5] and [1], 
but we shall also include some comments on the relationship between the isometries 
in the material space and in the spacetime. Most of the results in that section are 
well known and could be found in the above references, but we are still including 
them in order to set up the notation which will be followed in the remainder of the 
paper as well as for making the present paper more self-contained. Section 3 contains 
a digression on spherically symmetric spacetimes and the restrictions that such an 
assumption imposes on the physics in these spacetimes, which we then apply it to the 
case of elastic materials. In section 4 EFEs are obtained for the general case and some 
particular cases are commented upon. In section 5 we analyze in detail the case of 
shearfree solutions, paying special attention to the fulfilment of the dominant enrgy 
condition as well as to the necessary and sufficient conditions that must be satisfied 
for an equation of state to be admitted; finally we p resent a few selected examples; 
these include the analysis of the elasticity difference tensor for the nonstatic case, 
most along the lines followed in a previous paper by two of the authors [T3] . 

2 Relativistic elasticity revisited 

Let (M, g) be a spacetime, M then being a 4-dimensional Hausdorff, simply connected 
manifold of class at least, and g a Lorentz metric of signature (— , +, -|-, -|-). The 
material space X is a 3-dimensional manifold endowed with a Riemannian metric 7, 
the material metric; points in X can then be thought of as the particles of which 
the material is made of. Coordinates in M will be denoted as x"" for a = 0, 1,2,3, 
and coordinates in X as y^, A = 1,2,3. The material metric 7 is not a dynamical 
quantity of the theory, but it is frozen in the material, and it roughly describes 
distances between neighboring particles in the relaxed state of the material. 

The spacetime configuration of the material is said to be completely specified when- 
ever a submersion ip : M ^ X is given; if one chooses coordinate charts in M and X 
as above, the coordinate representative of is given by three fields 



and the physical laws describing the mechanical properties of the material can then 
be expressed in terms of a hyperbolic second order system of PDE. The differential 
map ip^: : TpM — T^^^p^X is then represented in the above charts by the rank 3 matrix 




A 



1,2,3 





A = 1,2,3, b = 0,1,2,3 



p 
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which is sometimes called relativistic deformation gradient. Since ip^, has maximal 
rank 3, its kernel is spanned at each point by a single timelike vector which we may 
take as normalized to unity, the resulting vector field, say u = u'^da, satisfies then 

= 0, u^Ua = -1, n° > 

the last condition stating that we choose it future oriented; u is called the velocity 
field of the matter, and in the above picture in which the points in X are material 
points, it turns out that the spacetime manifold M (or, to be more precise, an open 
submanifold of it) is then made up by the worldlines of the material particles, whose 
tangent vector is precisely u. 

The material space is said to be in a locally relaxed state at an event p G M if, at 
p, it holds kab = ('0*7)ab = hab where hab = gab + UaUb- Otherwise, it is said to 
be strained, and a measurement of the difference between kab and hab is the strain, 
whose definition varies in the literature; thus, it can be defined simply as Sab = 
— ^{kab — hab) = —^{kab — UaUb — gab)- Wc shall follow instead the convention in [3] 
and use 

Kab = kab - UaUb (l) 

Notice that K^u^ = u"", and therefore one of its eigenvalues is 1. Definitions using 
the logarithm of the above tensoi0 also appear in the literature as that allows simple 
interpretations of the associated algebraic invariants, see e.g. [T] and [S]. 

The strain tensor determines the elastic energy stored in an infinitesimal volume 
element of the material space (or energy per particle), hence that energy will be a 
scalar function of Kab- This function is called constitutive equation of the material, 
and its specification amounts to the specification of the material. We shall represent 
it as f = f (/i, /2, /a), where /i, I-z, h are any suitably chosen set of scalar invariant^ 
associated with and characterizing Kab completely. Following [3] we shall choose 

/i = ^(Tri^-4) 

l2 = ^ [TiK^ - (TtK)^] + 3 (2) 
/s = ^ (detA^ - 1) , 

Notice that for Kab = Qab (equivalently kab = hab) the strain tensor Sab is zero, that is: 
the induced metric on the rest frame of an observer moving with four-velocity u, h, 
coincides with the material metric 7 (its pull-back by ip) describing the relaxed state 
of the material; thus it makes sense to have zero elastic energy stored. It is immediate 
to check from the above expressions that in this case one has Ii = I2 = I3 = 0. 

The energy density p will then be the particle number density e times the constitutive 
equation, that is 

p = ev{h, h, h) = eo Vdet K v{h, h, h) (3) 

^ With one index raised, thus one has a Hnear operator which turns out to be positive and self- 
adjoint, its logarithm being then well-defined 

^ Recall that one of the eigenvalues is 1, therefore, there exist three other scalars (in particular 
they could be chosen as the remaining eigenvalues) characterizing if^ completely along with its 
eigenvectors. 



where eo is the particle number density as measured in the material space, or rather, 
with respect to the volume form associated with kab = {'4'*l)abi and e is that with 
respect to hab] see [TU] for a proof of the above equation. In some references (e.g. 
[3]), the names p and e are exchanged and the density measured w.r.t. kah = {'^*l)ah 
(eo in our notation) is then called "density of the relaxed material" (see the above 
comments on the meaning of 7), whereas that measured w.r.t. hab is referred to as 
the "density in the rest frame" . 

We next turn our attention towards the energy-momentum tensor of an elastic ma- 
terial. Before proceeding, it will be useful to recall that any symmetric, second order 
covariant tensor field may be decomposed with respect to a timelike unit vector field 
w, v"-Va = — 1 as follows: 

Tab = pVaVb + Phab + Pab + VaQb + QaVb (4) 
where hab = 9ab + VaVb, Pab = h"^hl{Tran - ^Phmn), Qa = -{TabV^ + P^'a), P = TabV^V^, 

p = ^h°'^Tab- From the definitions of these variables it readily follows 

habv'' = Q, Pabv' = g'''Pab = 0, and = 0. 

In the case that Tab represents the energy-momentum tensor of some material distri- 
bution, p, p, Pab, (f are respectively the energy density, isotropic pressure, anisotropic 
pressure tensor and heat flow that a family of observers moving with four-velocity v 
would measure at every point in the spacetime. 

In the case of elastic matter, it can be seen using the standard variational principle 
for the Lagrangian density A = \f--gp (see for instance [1] or [1] and the beginning 
of section 5 for further details) that the energy-momentum takes the form, when 
decomposed with respect to u, the velocity of the matter: 

Tab = PUaUb + phab + Pab (5) 

where all the definitions are the same as the ones given above substituting u for v, 
i.e.: hab = 9ab + UaUb, Pab = h"^h'l{Tr,^n " ^phmn) , P = TabWu^, p = Ih^^Tab and they 
satisfy habU^^ = 0, Pabu'^ = g'^^Pab = 0; thus in particular one gets g*^ = and the 
resulting tensor is of the diagonal Segre type {1, 111} or any of its degeneracies, u 
being its (unit) timelike eigenvector (see [TT]). 

This means that an orthonormal tetrad exists {ua-, Xa,yai Za} (with UaU'^ = — 1, x"'Xa = 
y'^Ua = z"'Za = +1 and the mixed products zero) with respect to which Tab may be 
written as 

Tab = pUaUb + PlXaXb + P2Z/aZ/fe + P^i^aZb, P = \^{Pl+ P2+ Ps) , 

hab = XaXb + VaVb + ZaZh, etC. (6) 

It is interesting to mention that the Dominant Energy Condition (DEC), see for 
instance [11], is fulfilled if and only if 

P>0, \pa\<P, a = 1,2,3. (7) 
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3 On symmetries and their consequences on physics 



Let (M, g) admit a Killing Vector (KV) C,, i.e.: in any coordinate chart x"", C^Qab = 0. 
It is then immediate to show that Rab = Gab = R%cd ~ 0' from 
EFEs it then follows that Tab = where T denotes the energy-momentum tensor 
describing the material in the spacetime. 

It is easy to show that, if iT is a non-degenerate unit (i.e.: v'^Va = e with e = ±1) 
eigenvector of Tab with corresponding eigenvalue A, then 

X = = 0. (8) 

We next include a short proof of this. 

(i) Taking the Lie derivative of TabV^, with respect to ^ and since we are 
assuming that Tabv'' = Xva one gets 

C^- {Tab V') = (A Va) = A C^- Va + VaC^- A. (9) 

On the other hand, 

C^^iTab v'')=Tab C^VK 

Equating (I9l) and ( fTOi) and then contracting with u", yields 
A Vb{C^ v^) = A v^'C^ Va + C^X. 

Therefore X = 0, since as ^ is a KV VaC^ v"" = v°'C^ Va- 

(ii) Substituting this result into ([9]) and using ffTOl) one obtains TabC^'^ = 
X C^a- Therefore, and the vector w"' = C^"" are eigenvectors of Tab 
associated with the same eigenvalue. Since this is non-degenerate, v and 
w have to be proportional, i.e. Cj^'"' = av"", for some real value a, however 

this a must be zero as v is unit and ^ is a KV, indeed 

= C^^Va = 2vaC^"' = 2av"'Va, hence a = 0. 

Under the hypothesis that is a Killing vector, and on account of the above consid- 
erations the following conditions hold in the case that Tab represents elastic matter 
and is therefore of the form ([5]) 



(10) 
(11) 



^^ab = ^ C^p = 0, C^a = 0, C^hab = 0, C^Pab = 0, = 0. 

The first two are just the specialization of the above comments to the case v = u and 
X = p. The vanishing of C^hab follows then from the vanishing of the Lie derivative of 
the metric and that of Ua', notice that C^h""^ = as well. Next, since p = ^h"'^Tab it also 
follows that its Lie derivative with respect to the KV vanishes as those of Tab and h""^ 
do, and the vanishing of Pab follows then immediately from the expression ([5]) and the 
vanishing of the Lie derivatives of all the other terms. Therefore we have shown that 



matter 4— velocity, pressure, density, anisotropic tensor all stay invariant along the 
Killing vectors of the space-time, together with the projection tensor hab = Qab + UaUb- 
It is interesting to notice that, for a general energy momentum tensor such as (jlj), 
and if one assumes that = 0, it also follows that = 0; but in this case that 
assumption has to be made, as u is no longer an eigenvector of the energy-momentum 
tensor. 

Similar conclusions (although not the same) can be drawn when ^ is a proper ho- 
mothetic vector field; i.e.: gab = '^kgab, with /c 7^ 0; in which case one also has 
-^1* Tab = but then, for instance, Ua = —kua and also p = kp, etc. 

In this paper we shall be concerned with the case of elastic materials in spheri- 
cally symmetric spacetimes. We next explore the consequences that the existence of 
symmetries has on the material content of a spacetime. Most of the developments 
following are well known although disperse in the literature, we collect them here for 
the sake of completeness. 

As it is well known, for a spherically symmetric spacetime, coordinates x° = t, r, 6, (f) 
exist (and are non-unique) such that the line element can be written as 

ds^ = -a(r, t)dt'^ + b{r, t)dr^ + r'^dO'^ + sin^ ed<f (12) 

with a and h positive and independent of 9 and (j). This metric possesses three 
Killing vectors, namely ,^1 = — cos0 de + cot 6^ sin 9^, ,^2 = and ,^3 = 
— sin (j) dg — cot 6 cos (p which generate the 3-dimensional Lie algebra so(3). 

To start with, we show that any timelike vector field v that remains invariant along 
the three Killing vectors is necessarily of the form 

V = w*(t, r) dt + f^(t, r) dr- 

Using jC^^v'^ = for a = 0, 1,2,3 we conclude that all the components v"" are inde- 
pendent of (p. Then, the expression 

for a = 0, 1 gives that and are also independent of 6 and for a = 2, 3 yields 
v'^ = = 0. 

It should be noticed that it always exists a coordinate transformation taking t, r into 
t', r' such that v = {t', r') df and the metric (fT2l) reads then 

ds^ = -air, t)dt^ + 6(r, t)dr'^ + ^^(r, t) {dO^ + sin^ ed(j)'^) (13) 

where primes have been dropped for convenience. These new coordinates are the so 
called comoving coordinates; one can see this either by direct computation, showing 
that one such coordinate change is always possible, or else, by showing first that 
any vector field such as v above is always hypersurface orthogonal, that is: uiab = 
V[a;b] + V[aVb] = 0, then it follows that Va oc dat' for some function t', choosing this 
function as the new time coordinate, one can readily show the above result. 

Next, for any symmetric, second order tensor Pab, which is traceless, invariant under 
the three KVs above, and orthogonal to a vector such as v above (i.e.: timelike 
spherically symmetric), that is: 



fi 



(i) Pabv' = 
(ii) g^'Pab = 

(iii) %P„, = 0, for A = l,2,3 

it follows that Pab is proportional to the shear tensor of v whenever the latter is 
non-zero, namely: 

Pab OC aab, (Tab = V(^a;b) + V(^aVb) " ^Ohab (14) 

where ^ = w!^ is the expansion, Va = Va-bv'^ is the acceleration and round brackets 
denote symmetrization as usual. 

This can be proven easily by making use of the comoving coordinate system referred 
to above, and imposing the various conditions (i - iii); also, a more general proof is 
possible in the context of warped spacetimes (of which the spherically symmetric ones 
are special instances), see [T2] . 

In the comoving coordinate system above, one can see by direct computation that 
the shear aab is 

/ l{htY-2Yth) lYihtY -2Yth) lYihtY -2Yth) . \ , , 

and therefore this field is shearfree if and only if 

h{r,t) = F\r)Y\r,t), 

in which case it is always possible, by means of an obvious redefinition of the coordi- 
nate r, bring the metric to the form 

ds^ = -a(r, t)dt^ + Y\r, t) {dr^ + dO^ + sin^ Odcf)'^) . (16) 

For the class of spacetimes we shall be interested in, namely elastic, spherically sym- 
metric, all the above apply for u (the velocity of matter), as it is indeed the unit 
timelike eigenvector of Tab given by ([5]), and Pab the anisotropic pressure tensor since, 
as discussed previously, it is invariant under the KVs the spacetime possesses and is 
also traceless and orthogonal to u, therefore we have that, whenever the shear of u is 
non-zero 

Pab = 2X(Tab, (Tab = U(a;b) + ^^(0^6) " ^^^afe (17) 

where 9 = u""^ and iia = Ua-bu'^, and A = A(t, r) is some function, therefore, for the 
generic (non shearfree) case, it is always possible to treat, at least formallj{§, the 
elastic material as a viscous fluid with zero heat flow. This interpretation would 
indeed break down in the case in which u is shearfree. 



^ There is a further requirement for a fully physically meaningful interpretation as a viscous fluid, 
namely that A < in which case the kinematical viscosity would be i] = —A > 0. 
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4 Elasticity in spherical symmetry 



Let us now consider in more detail the problem of elasticity in a spherically symmetric 
spacetime {M,g) with associated material space (X, 7). 

The results given in this section generalize those in [3] in the sense that here we 
consider a non flat material metric 7, while, when referring to quantities and results 
in [3], we shall use non-barred quantities (hence the bars on the spacetime metric and 
the material metric in our notation). 

Recalling the notation and results in section 2, we shall demand that the submersion 
ip : M — > X preserves the KVs, that is: V'*(Ca) = Va are also KVs on X. 

This implies that the metric 7 is also spherically symmetric and therefore coordinates 
= (2/5 ^5 0) exist with y = y{t, r), 9 = 9 and (f) = 4>, and are such that tJa = C,a are 
KVs of the metric 7. Thus, the line elements of g and 7 may be written as: 

df = -a{t, r)dt^ + h{t, r)dr^ + r^d9^ + r^sin'^9d(f)'^ (18) 



rfS^ = f{y){dy^ + y^d9^ + yhin^9d^^), (19) 

Notice that this last expression is completely general, as any 3-dimensional spherically 
symmetric metric is necessarily conformally fiat, as it is immediate to show. 

The results in |3] correspond to f{y) = 1, and the relation between 7 and the fiat 
material metric 7 used in |3] is given by 

lAB = f{y)lAB. (20) 



Next, attention should be payed to the canonical definition of the energy-momentum 
tensor used by [5]: 



rpa 

-'-b 



dA 



9 \dy^ 



(21) 



As shown in [T3], this canonical definition of the energy- momentum tensor coincides 
with the symmetric definition of the energy-momentum tensor, used by other authors, 
up to a sign, which is a particular case of the general Belinfante-Rosenfeld theorem 



Denoting by k the pull-back by ip of the material metric 7, that is: k = ip* one 
has: 



C„.B 



cb 



ny)r[fsX + yy'{5X + 5^) + yX^l + yX^l + y's^nV,sf 



cJ' 



or 



( -p{y){yVd) -P{y){yy'/-a) \ 

P{y){yy'/b) P{y){y'Vb) o o 

P{y)yVr^ 

V P{y)yVr' j 



(22) 



where a dot indicates a derivative with respect to t and a prime a derivative with 
respect to r. 

The velocity field of the matter, defined by the conditions u°'ya = 0, gabU°'u'^ = —1 
and > 0, can be expressed as 



where 



r 



1,-^,0,0 

y 



a \y' 



(23) 



(24) 



Therefore the projection tensor is 

/ l-f2 



\ 



-T\by/ay') \ 

F(W) l+T'ib/-a)iy/yr 

10 

1 / 



(25) 



We will use an orthonormal tetrad and write the metric as Qa 
yaVb + ZaZb, such that: 



X" 



r y T 

Va y' Va 

-^Ir L 

a y Vb 



0,0, -,0 
r 



Ua 



-UaUb + XaXb + 



ar,-4=^f,o,o 
vay 



^v^r, v^r,o,o 



0,0,0, 



1 



r sin 9 



ya = {0,0, r,0) 
Za = (0, 0, 0, r sin i 



where T is the auxiliary quantity given in (l24ll . Here, is the matter velocity and 
x", y"" and are spacelike eigenvectors of the pulled-back material metric k^. From 
our developments in section 2, it is immediate to see that the pressure tensor has the 
same eigenvectors as kab and can be written, for the space-time under consideration 
as Pab = piXaXb + P2{yayb + ZaZb) . Therefore, ^ yields 

fab = pUaUb + plXaXb + p2 {VaVb + ZaZb) , (26) 

where p is the energy density, pi, the radial pressure and p2, the tangential pressure. 

The results in |3] can be easily recovered by setting f{y) = l above. 

Now, much clarity is gained by making use of the comoving coordinates adapted to 
u, the timelike eigenvector of the energy-momentum tensor, which were introduced 
in the above section. The form of the metric is given by 



u, being then 



-a(r. 



, t)dt^ + b{r, t)dr^ + Y\r, t) {dO^ + sin^ ed(f)^) 



^'^=(^-^,0,0,0^, = (-v^, 0, 0, 0) 



(27) 



f28) 



hence, we have for the material space (M, 7) that coordinates = (y, 9, 0) exist 
with y = y{r), 9 = 9 and = 0, as follows from the condition y^u"" = and the 
requirement that tp^:{C,A) = Va are KVs of the metric 7. 

Further, and since the line element of the material space is 

with y = y{r), no generality is lost by setting y = r, as this amounts to a redefinition 
of the r coordinate in spacetime, and leaves unchanged the form of the metric (127|1 
as well as that of the velocity field of the matter (128|) . We shall do that in the sequel. 



Thus, the pulled-back material metric k 



IS 



where a prime indicates a derivative with respect to r, which upon setting y = r as 
discussed above it simplifies further to: 



/ \ 

/2(r)(l/6) 

Piry/Y"^ 

\ p{ry/Y^ J 



(29) 



The operator = g'^^kcb — u^-Ub, introduced in section 2 and used to measure the 
state of strain of the material has one eigenvalue equal to 1 (corresponding to the 
eigenvector u), while the other eigenvalues are 



f{y) = fir) 



/2 



n = fiv) %- 



(30) 



and s has algebraic multiplicity two. 



The three invariants /i, /2, 13 of K introduced in ([2]) have the following expressions 



Ji = ^(Tri?-4) =^(r/ + 2s-3) 



h 



liK'^ - (Tr/i ) 



4 

\ (deti^ -l) = \ ivf - 1) 



3 = {s^ + 2f]s + f] + 2s) -3 



(31) 



In [3], the energy-momentum tensor was calculated from these invariants for a flat 
material metric. A similar calculation shows that, for the non-fiat material metrics 
under consideration, the same expression holds so that 



dp 



- dp dp 



= pSl--J^ detK hi + TrA' - ^ ] kl - ^ k"^ k^ 



dh dh 



9P La 



dh 



(32) 



in 



Therefore, the nonzero components are 



- n 



,'2 



fl = fl = p- 



Y2 



9p 



f \y) 



dp 



y2 dh 



7/2 7/'2 



where 



The rest frame energy per unit volum^, p, is defined by 

p = ev= eos ^v{s,f]), 



(33) 



(34) 



(35) 



where, as discussed in section 2, v = v{Ii, h, h) = v{s, ff) represents the constitutive 
equation, eo, the density of the relaxed material (density w.r.t the pulled-back material 
metric k) and 



e = eo V detA' = eo s y^, 
the density calculated in the rest frame (that is, w.r.t. h). 
Then, using (13 ip . one can prove the following relations: 
dp 1 



(36) 



ds P 



dh 



dp_ 
F2 dh 



Y2 



.12 



(37) 
(38) 



Alternatively, one can express the components of the energy-momentum tensor in 
terms of the eigenvalues s and fj by substituting the last results in fl55]) : 



To" = e-7;, 

dv 
ds' 



(39) 



-'2 



-e s 



The Einstein field equations Gl = SttT^" can be written as follows: 



Y 



Y b 2Y" Y 



12 



Y' h' 



YH Yah Yb Y% Y Y^ 



(40) 



2Y' --Y- Iy' = 0, 
a b 



In [3] the quantities p and e are e and p, respectively. 



(41) 



1 1 



y2 Y a Y' a' Y'^ 2Y 1 dv 

+ ^r^ + ^TT + 7>;T-^-^ = -e2r/ — Stt, (42) 



Ya? Yah Y% Ya ' 
Gl = Svrfl: 

1 Fa 1 FS 1 a'2 1 F'a' 1 Y" 1 F'fr 



2ya2 2Fa6 4:a?h 2Yab Aab^ Yb 2Yb^ 
la" li 1 M> IP f _ 

2dS~2ub^ Am ^ 4^ ~Yi~ 

__dv 
-es — STT. 

OS 



(43) 



It is interesting to express the contracted Bianchi identities for in terms of v and 
its derivatives w.r.t the quantities f] and s. Thus, from T^.^ = one has: 

fl, + 9„(ln v/^)f," - Tlf: = (44) 

and specifying this equation to 6 = 0, 1 one gets respectively (for non- stationary 
solutions): 

dtiev) + ^ ^-ev + -ev-j + ^{2ev + 2e.--) = 0, (45) 



/ dv\ 1 a' dv fa' Y'\ dvY' 

The remaining equations for 6 = 2,3 which can be obtained from (jHj) are identically 
satisfied. 



Equation (1451) for non-stationary solutions, implies readily 



.-^.„W, (47) 



which can then be substituted into (1461) to get a slightly simplified equation. 

From this point onwards, we shall drop the bars, as no confusion may arise with the 
results in [3]. 



5 Shearfree solutions. Examples 

In this section we shall consider in detail the case of spacetimes with a material content 
that may be represented by some elastic material such that the velocity of the matter 
is shearfree, in which case coordinates exist such that the metric can be written in 
the form (lT6l) . For this case, the interpretation as a viscous fluid with kinematical 
viscosity is not possible, and therefore the anisotropy in the pressures must be a 
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consequence of the elastic properties of the material. The study of solutions with 
non-vanishing shear tensor and their possible interpretations as viscous fluids, will be 
carried out elsewhere as this would render the present paper too lengthy. 

We will study separately the cases of static and non-static solutions, presenting ex- 
amples of each instance which are regular at the origin, posses an equation of state 
and satisfy the dominant energy condition (at least in some open submanifold of the 
spacetime). 

Consider the metric (fT6|) which we rewrite here for convenience: 

ds^ = -a{r, t)df + Y\r, t) [dr^ + dO^ + sin^ ed(j)^) (48) 
From the field equations it follows that G*^ = which in turn implies that 

a = L{t)y (49) 
whenever Y ^ 0, L{t) being a function of time. 

If y = then G^^. = is identically satisfied, and from (!30|) follows that s and rj, and 
therefore v{ri, s) are functions of r alone, then (145!) implies that e = e(r); further from 
the field equation G\. = — 87re2?7|^ it follows that G\. can only depend on r as well, 
which in turn implies that a{t,r) = ao{t)ai{r), the solution being then static, as a 
trivial redefinition of the coordinate t coordinate shows. 

It is interesting now to see that in the shear-free case, if one sets either r) or s equal to 
1, so that matter is strained in tangential directions (but not in the radial direction 
?7 = 1), or it is strained only in the radial direction [s = 1), from the definition of these 
quantities it follows that Y = Y{r), and according to the statements in the above 
paragraph, it follows that the solution must be static, and therefore the results in [8] 
apply. Thus, we have proven that: if the velocity field of the matter is shear-free and 
the matter is stressed either in the radial direction only or in the tangential directions 
only, the spacetime is necessarily static. 



5.1 Static shearfree solutions 

In the static, shear-free case (metric (HHj) with no dependence on t), the field equations 
yield 

Y" y2 
dv a'Y' Y'^ 1 
dv Y" 1 a" _ 1 a'2 ^ y/2 

ds^"" ~ Y'^ ^ 2 aY' Aa'Y^ ' (52) 

solving (150!) for e and substituting it in fl5T|) and fl52l) one gets two equations which 
depend only on r and elementary considerations show that for given a(r) and Y{r), 



functions y{r), f{y) and v can be found so that the two equations are satisfied. It 
remains to be seen, though, that the DEC are satisfied and therefore the solution is 
physically acceptable. 

The following simple example shows that solutions with these characteristics do indeed 
exist. 

Example 1 

Consider the line element 

ds^ = ~Y-\r)dt'' + y'(r) {dr^ + de"" + sin' Odcj)^) . (53) 
A direct calculation yields 

The dominant energy condition (JTj) implies: 

Y" y2 ^ 

8^P = V 7^"!^-°' ^^^^ 
87r(p-|>i) = 2^>0, (56) 

/Y" I \ 

8^(P + ^'i) = 2(^-yi-^j >0 (57) 
8^(p-P2) = 2^-2^-i^>0, (58) 

87r(p + p2)=2^-i^>0, (59) 
where we put p = ev, pi = —2er]^ and p2 = ^^-^f^- 

Now, it is immediate to see that the above conditions are all satisfied if and only if 
(l57|l is, which in turn can be written as 

y-\-y-y^-^)>^ ^ (lnF)"-l>0, (60) 

which is equivalent to 

F = exp (rV2) /2(r) such that (In/)" > 0. (61) 

Take, for instance, 

Y = exp(5/2r2), (62) 
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one then has 

p = ev = -^e-^"'(25r2 + 9), 

OTT 

Pi = -2er/^ = -^e-^'-'(25r2 + 1), = -es^ = ^25r'e~''-' (63) 

or] Sn OS Svr 

which is obviously well behaved: satisfies the dominant energy condition and is non- 
singular at the origin. Notice that the radial pressure is negative (compressed mate- 
rial) and the tangential pressures are zero at the centre, as one would expect. 

The field equations in this case read: 

ev = —e-^''\25r^ + 9) (64) 
Svr 

- 2er]^ = -— e-^^'(25r2 + 1) (65) 
orj Sir 

_ = Le~^r^25r^ (66) 
ds 871 ^ ' 

and one has that 

r^ = f^r)e-''\ s = p{r)e-^'\ (67) 

Now, dividing fl65l) and fl66|) through by flM|) . and setting £" = ln?7, E = In s, one gets 
dhiv 125r^-hl dlnv 25r^ 



dE 2 25r2 + 9' 25r2 + 9' 



(6J 



From the expressions for rj and s one has that 

E = 21n/(r) -5r2, E = E + 21nr, (69) 

hence one can express r as a function of and E as a function of E as well, thus 

d In V (9E dh\v 
dE ~ dE 9E 

from where it follows that 

9r 3 1 2 

M = -r-1007' i? = -3ln(75r^ + l), 

that is 



e-i^-1). (70) 



75 

Plugging the expression of E in terms of r into (l69l) one gets that 



(75r2 + 1)3 
1.1^ 



whence expressions for 77, s and e = eQSy/rj can be easily derived. 

Next, from (1701) and the first equation in ( l68i) . one can easily find an expression for 
the equation of state, namely: 



where -F(S) must satisfy 



e2^ 
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(72) 



ainF 25r^ 



dE 25r2 + 9' 

where r in the right hand side of the equation has to be expressed in terms of S. 
From the second equation in fl69l) it follows that r must be the only real solution of 



which has a rather complicated form. In any case, one gets 

F(S(r)) = (25r2 + 9)"^ (TSr^ + l)~ 



_L_ 

39 



and thus we have proven that a solution exists, which is regular at the origin r = 0, 
satisfies the dominant energy condition and possesses an equation of state which can 
be given in a closed form. 



5.2 Non-static shearfree solutions 

Assume now that Y ^ 0, so that a{t,r) takes the form (H9!) . substituting this into 
(I48p and redefining the coordinate t so as to absorb the arbitrary function L(t) one 
has 

ds^ = -Ldf + Y^ (r, t) [dr^ + dO^ + sin^ Odcf)'') (73) 

From G*^ = it follows now that F(r, t) = A(t)B{r), which substituted above, yields, 
after a trivial redefinition of the coordinate t: 

ds^ = -de + A^{t)B^{r) [dr^ + dO^ + sin^ ed(P'') (74) 

A direct computation of the EFEs for the above metric give gives 

_ 2B"B + ?>A^B^ + 52 
Svret; = , (75) 

dv 2AB*A-B'^ + A^B* + B^ , , 

- ^^-^^Q-^ = J^. ' (76) 

dv 2AB^A- B"B + A^B^ + B'^ , , 

- ^^^^^ = • 



where the energy density, radial and tangential pressures are 

^ dv dv 
p = ev, pi = -2e?7— , p2 = -es—. 

OT] OS 

On the other hand, the dominant energy condition, p>0,p±pi>0 and p ± p2 > 
implies 



- - 2B"B + SA^B^ + B^) > 0, (78) 

- B'^ + B"B - A^B^ + AB^A > 0, (79) 
- [AB^A - BB" + 2A^B^ + 5^) > 0, (80) 

BB" - 2A^B^ - B^ + 2AB^A > 0, (81) 

- (25'2 - 3B"B + AA^B^ + B^ + 2AB^A) > 0. (82) 

We next address the question of the existence of an equation of state v = v{'r],s) 
where in this case, t] and s are given by (see (130|1 ) 



' A2(t)52(r)' A^{t)B^{r)' 

Dividing fl76p and (1771) by (175]) . and defining as before = ln?7, E = Ins, we get 
dlnv 1 2AB^A - B'^ + A^B^ + B^ 



dE 2 _ 2B"B + 3A^B'^ + B^ 



Ve, (84) 



dlnv _ 2AB^A - B"B + A^B^ + B'^ _ ^ 
OT. ~ 5'2 - 2B"B + 3i254 + 52 " ^' 

In order for an equation of state v = v{r], s) (or equivalently v = v{E, S)) to exist, it 
must be that 

dWE ^ dEdT. ^ 'OT^^'dE' ^ ' 

Notice that 

Now, from the expression (l83l) and the corresponding one for E and S, it follows that 

r = e^(^-^), A{t) = 4e-^^; (88) 

B 
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and differentiating tliem witli respect to E and S and applying tlie cliain rule, we get 

(89) 



dr 
dE 



1 

2 ' 



dr 

as 



1 

-r, 
2 ' 



dt_ 
dE 



1 

2A 



/ 



AB dt _ 1 ffVAB 
"T' ^"21'' [bJ ~T' 



(90) 



Substituting the above expressions into (1871) . equation (l86l) reads, after some manip- 
ulations 



i (y " f") + + ^ -^-(^^ + Vs) + j dtVs = 0. 



(91) 



/' 

Solving for — we get, after some algebra. 



2,rB 



2f 
f 



-K 


{S + ?>B^A^)Attt - 6B^AA^ + 3B*A-^A^ 


+ MA-^AA 




{S + W^A^)Aut - QB^AA^ + 3BM-ii3 


+ SA-^AA 



(92) 



where, K, S and M are functions of r alone given by: 

K = 2B- 3rB', S = B'"^ - 2BB" + fi^ 
M = -QrB^Brrr + 2B{9rB' + B)fi" + (2fifi' - 9rB'^ - 3rB^)B' - AB'\ 

Since the left hand side of (192|) depends only on r this implies that the time derivative 
of the right hand side must vanish. A careful but otherwise trivial analysis, reveals 
that there are only three possibilities, assuming the metric is n on-static, namely 



1. M + KS = 0, in which case B{r) is determined by the resulting ordinary third 
order differential equation. A{t) is in principle arbitrary, and /(r) is fixed by 
(p2l) once the solution for B{r) to the equation M + KS = is give n. We have 
not been able to find an integral for B{r) in closed form, but in this case one 
has for /(r): 

f _ B' _ 2 
~~B^ 



/ 



3rB 



2. A = 0, which in turn implies A = t without loss of generality, since the two 
constants of integration may be absorbed by suitable redefinitions of t and B. 
In this case, B{r) is free, constrained only by the requirements imposed by t he 
DEC, and once it is chosen, /(r) is determined through fl92|) . which implies as 
in the previous case 

f _ B' 2 
J ~B^~ 3^' 



1R 



3. In this case, both A(t) and B{r) are determined as the solutions of the following 
two third order differential equations: 

-2BB" + B'^ + B^ - kB"^ = 0, k = constant, 
Since A{t) and B{r) are fixed, so is f{r), and fl92l) implies in this case 



,f' _ M + {k-3q)KB^ 
J ~ ■ 

As in the first case above, we have not been able to find integrals for A{t) or 
B{r) in closed form. 

Example 2 

Let us next investigate in some detail the second case above, that is A = t. Substi- 
tuting this into the EFEs we get 

1 f2B" 5'2 1 \ _ 1 1 \ _ 1 fB" B'^ 

\w ~^'b^~^)' P'-s^ym^B^^y' ^'~8^[b^~^~^i -^^^ 

On the other hand, the DEC (178 i) -( l82l) are all satisfied if and only if the following 
three inequalities hold: 

3BB" - 2B'^ - B^- AB^ > 0, BB" - B'^ - B*>0, BB" - B^ - 2B^ > (94) 
which, upon setting B = are equivalent to 

3b" + b'^ - 1 - 4e^' > 0, 6"-e2^>0, b" + b'^ - 1 - 2e^' > 0. (95) 

It is easy to see that these conditions can be satisfied, at least for certain fanges of 
the radial coordinate r G [0,R), for suitably chosen functions b{r), such as 



b{r) = 2 



— h sinh^ f 

3 \ c 



Inc -v^ B(r) = -^In 
^ ' 9c 



-1 + 3cosh2 



2 / r - ro 



(96) 



The form of /(r) can be given explicitly up to a quadrature. 

While we do not claim that it has any particular significance, it provides a relatively 
simple instance of solution with the desired properties. 

Example 3 

Another simple example with similar characteristics is provided by the following choice 
of 5fr): 



%) = ^lnQ + r^) ^ B{r) = y^{2 + 3r')\ 

1Q 



(97) 



in which case /(r) can be integrated out yielding: 



/(r) = exp 



-15 - 9r2 + 1(2 + 3r2)V4 + 6r2tanh"^ {vm^) 
V6 + 9r2 (2 + 3r2) 



Similar remarks to the ones in the previous case regarding its physical significance, 
apply also here. 



6 The elasticity difference tensor for non-static solutions 

Here we obtain the elasticity difference tensor, defined in [1], for non-static spherically 
symmetric spacetimes and analyze this tensor following the procedure developed in 
|13j . where the static, spherically spacetime case was presented as an example. 



This third order tensor, symmetric on the two covariant indices, is completely flow-line 
orthogonal and is related with the (pulled back) material metric according to 

S\c = \k-''"'{Dtk„,, + D.kmb - (98) 

Here k'"'"^ is such that k~°'"^kjnb = and D represents the spatially projected con- 
nection obtained from the spacetime connection V associated with g by 

Djt:: = hthl..hl..v,ti-:, m 

with the property Dah{,c = 0. 

The non zero components of S^^ for non-static, spherically symmetric spacetimes, 
using the space-time metric (l27|l and the pulled-back material metric (129|) can be 
written as: 

f b' 



rr 



C<P 
or 
or 



/ 2b 

f 1 y 

f^2 ^ Y'Y 

_i \- 

f r b 

r^f sin^ 9 r sin^ 9 Y'Y sin^ 9 



f 1 



In this case, the pulled-back material metric kab is 

kab = n\ XaXb + nj {yaVh + ZaZ^). (100) 

Here, x, z are eigenvectors of kl with eigenvalues n\ and n\ = which depend on 
t and r according to 

nl = fl nl = nl = f^^. (101) 



90 



The elasticity difference tensor can be decomposed along the directions determined 
by the eigenvectors of as follows 

= + M,,y'^ + M,,z\ (102) 

12 3 

where M, i = 1,2,3 are second order, symmetric tensors (see It should be 

i 

noticed that the eigenvectors 

Here we determine the eigenvectors and eigenvalues of these tensors, complementing 
the results obtained in [13] for the static case, the result being summarized in Tables 
1,2 and 3: 



Table 1 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


X 

y 

z 


u - ^" 

Vbf 2bVb 
Y' r^f'Vb rVb 

/^3 = 



Table 2 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


X + y 

x-y 

z 


/' 1 y' 

^4 " TTf + ~ FT! 
/i5 = -/^4 
yue = 



Table 3 - Eigenvectors and eigenvalues for M 



Eigenvectors 


Eigenvalues 


X + z 
X — z 

y 


/i7 = /i4 
12s = -/i4 
/ig = 



Therefore, the canonical forms for the three tensors M are: 

i 

Mbc = fJ'lXbXc + f22{ybyc + ZbZc) 

1 

Mbc = 2fi4{xbyc + ybXc) (103) 

Mbc = 2n4{XbZc + ZbXc). 
3 

Although the eigenvalues are different from the ones obtained in the static case, the 
eigenvectors of the above tensors are the same for the static and non-static case. 

7 Conclusions 

In this paper we have considered spherically symmetric spacetimes with elastic mate- 
rial content . We started considering the symmetries these sapcetimes posses in order 
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to exploit their consequences on physics for elastic spacetime configurations. By doing 
this, we have generalized previous work done by [1] and [3] for non-static spherically 
symmetric configurations, where only fiat material metrics were considered. In fact, 
we have shown that all material metrics compatible with a given spacetime are con- 
formally related and, moreover, are conformally fiat. Next we have used comoving 
coordinates to relate the EFEs with quantities characterizing elasticity properties 
(constitutive equation, material and energy density, eigenvalues of the pulled back 
material metric) as well as the conformal factor referred to above. The case in which 
the velocity of the matter is shearfree has been considered in detail, giving the neces- 
sary and sufficient condition for a constitutive equation be admitted; further, we have 
provided three examples of static and non-static shearfree solutions. For non-static 
spherically symmetric space-times, the elasticity difference tensor has been studied, 
thus extending some previous work for the static case. 
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